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Using the finite-size scaling, we have investigated the percolation phase transitions of evolving 
random networks under a generalized Achlioptas process (GAP). During this GAP, the edge with 
minimum product of two connecting cluster sizes is taken with a probability p from two randomly 
chosen edges. This model becomes the Erdés-Rényi network at p = 0.5 and the random network 
under the Achlioptas process at p = 1. Using both the fixed point of s2/sı and the straight line 
of ln sı, where sı and s2 are the reduced sizes of the largest and the second largest cluster, we 
demonstrate that the phase transitions of this model are continuous for 0.5 < p < 1. From the 
slopes of Ins; and In(s2/si)’ at the critical point we get the critical exponents 6 and v, which 
depend on p. Therefore the universality class of this model should be characterized by p also. 
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The modern theory of complex networks has two connecting cluster sizes is taken with a probability 


opened new perspectives in the study of complex sys- 
tems in nature, social and economic systems, technical 
infrastructures and many other fields. The macroscopic 
properties of complex networks emerge from the inter- 
actions among individual constituents. The percolation 
phase transition of complex networks is an interesting 
macroscopic property and can be studied with the per- 
colation theory in statistical physics It has been 
pointed out that the random network undergoes a con- 
tinuous percolation phase transition during the random 
process |5]. The critical phenomena of complex networks 
are reviewed in Ref. |6]. 


However, Achlioptas et al.{7] reported that the perco- 
lation phase transition of random networks becomes dis- 
continuous under the Achlioptas process (AP). During 
this AP, two unoccupied edges are chosen randomly and 
the edge with minimum product of the connecting cluster 
sizes is taken as the next occupied bond. This discontin- 
uous phase transition in networks is called the explosive 
percolation i. The AP was applied later to other net- 
works and some similar rules were also introduced 
[H15]. Using the finite-size scaling, the explosive per- 
colation phase transition was also investigated 
These works support that the explosive percolation is 
discontinuous. Recently, Costa et alg showed that the 
explosive percolation transition is actually a continuous 
phase transition with a uniquely small critical exponent. 
Riordan and Warnke prove mathematically that all 
Achlioptas processes have continuous phase transitions. 


In this letter we investigate the percolation phase tran- 
sitions of random networks under a generalized Achliop- 
tas process (GAP). In the GAP, two unoccupied edges are 
chosen randomly and the edge with minimum product of 


p 22. At p = 0.5, this model becomes the Erdös-Rényi 
(ER) network [25]. At p = 1, the GAP is equivalent 
to the Achlioptas process which suppress the appearance 
of larger clusters. With the increase of p from 0.5 to 1, 
this suppression is switched on gradually. After inves- 
tigating the percolation phase transitions at different p, 
we can get more understanding about these phase tran- 
sitions and make a more reliable conclusion about the 
character of the percolation phase transition at p = 1. 

In our Monte carlo simulations, we begin with N iso- 
lated nodes and then connect them with the edges added 
through the GAP. The network obtained is character- 
ized by N, the number of edges N, and the probabil- 
ity parameter p of GAP. For p < 1, we take N = 
(4, 8, 16, 32, 64, 128) x 10* and run 500, 000 steps in 
each simulation. At p = 1, we choose much larger system 
sizes N = (512, 1024, 2048) x 10* to reach the asymp- 
totic region of finite-size effects and run 100, 000 steps in 
each simulation. 

For the cluster ranked R and with size Sp(N,, N; p), 
we define its reduced size as 


sp(r, N; p) = Sr(N,, N; p)/N , (1) 


where r = N,/N is the reduced edge number. The 
reduced size of the largest cluster si(r,N;p) at N = 
1.28 x 10° is shown in Fig. [I] sı becomes finite for r > re 
and there is a percolation phase transition. The tran- 
sition point re increases with p and the corresponding 
phase transition becomes sharper. 

If the percolation phase transition is continuous, the 
reduced sizes sr(r, N;p) should follow a finite-size scal- 
ing form 


sr(r, N; p) = NP 3RN"; p), (2) 


Figure 1: 
10°. 


Reduced size of the largest cluster at N = 1.28 x 


where t = (r—re)/re and v is the critical exponent of the 
correlation length € = ĉo |t| ”. The scaling form in Eq. 
(2) is valid in the asymptotic critical region with |t| « 1 
and N œ 1. Outside this region, additional correction 
terms should be taken into account. 

From Eq. @), we can obtain 


82/81 = 32(tN/”; p)/ (tN; p) =U(tN;p) (3) 
and 
Ins1(r, N;p) = —B/v nN + n31 (tN; p). (4) 


At a critical point, s2/s1ı|r=r. = U (0; p) is a fixed point 
versus N and lIn sı (re, N; p) = —8/vln N + ln 51 (0; p) is 
a straight line versus ln N. Using these properties, the 
critical point of complex network can be determined both 
from the fixed point of s2/sı and the straight line of 
ln sı versus In N. The critical exponent ratio 8/v can be 
obtained from the slope of ln sı at re. 

For p = 0.5, our model is equivalent to the ER model 
fos}. Its continuous percolation phase transition can be 
studied with s2/s, and ln s1. In Fig.BYA), s2/s1 is shown 
as a function of r at different N. The different curves of 
89/81 for different N have a fixed point at rt = 0.5000(4). 
In Fig. B{B), Ins; is shown versus In N for different r. 
The curvature of ln sı is negative at r = 0.4996 and be- 
comes positive at r = 0.5004. At r = 0.5, Ins; is a 
straight line with slope —0.33(1). We get r? = 0.5000(4) 
in consistence with rl. Our results of re and 8/v agree 
with the exact results re = 0.5 and 8/v = 1/3 of the ER 
model (25). 

In Fig. B| we show s2/s; and ln sı of random network 
under the GAP at p = 0.8. From the fixed point of 
82/81, we get the critical point ri = 0.7252+0.0001. The 
curvature of ln sı is negative at r = 0.7250 and becomes 
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Figure 2: (A) s2/s1 has a fixed point at r4 = 0.5000 + 0.0004; 
(B) Log-log plot of the reduced size sı versus N. Around 
r2 = 0.5000, the curvature of In sı changes its sign. 


positive at r = 0.7254. So rł = 0.7252 + 0.0002, which is 
consistent with rl. 

For p = 1.0, our model becomes the PR model of 
Ref. (7. 82/8, at different N are shown in Fig. AA) and 
there is a fixed-point at r? = 0.88844(2). The curvature 
of In sı is negative at r = 0.88840 and becomes positive at 
r = 0.88850. We obtain r2? = 0.88845(5) consistent with 
ri. From the slope of ln sı at r2, we get 8/v = 0.04(1). 
Our re and 8/v at p = 1 agree well with Ref. (26). 

To determine the critical exponent v of correlation 
length, we study the derivative 


(so/ si)! = ŽEN L ety” sp) (5) 


At a critical point, we have 
In(s2/s1)'|rar, = (1/v) n N + InU'(0;p) —Inre, (6) 


which is a straight line with slope 1/v. From our Monte 
Carlo data, we can calculate 1/v at different p, which are 
summarized in Table I. 
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Figure 3: (A)s2/s1 has a fixed point at r4 = 0.7252 + 0.0001; 
(B)Log-log plot of the reduced size sı versus N at different r. 
Around r2 = 0.7252, the curvature of In sı changes its sign. 


Table I: Critical reduced edge number re, critical exponent 
ratio G/v and inverse of critical exponent 1/v. We obtain r2 
from s2/s1 and r2 from ln sı. 


Using re and 1/v given in Table I , we can define a 
finite-size scaling variable tN‘/” at each p. With this 
scaling variable, our Monte Carlo simulation data of 
s2/sı(r, N; p) at a given p and different N collapse into a 
finite-size scaling function U(tN'/”; p), which are shown 
in Fig. [5] for p = 0.5 and p = 1. The p-dependence of 
U(0;p) is demonstrated in Fig. [6] 
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Figure 4: (A) s2/sı has a fixed point at ri = 0.88844 + 
0.00002. (B) Log-log plot of the reduced size sı versus N. At 
r2 = 0.88845, the slope of ln sı versus In N gives B/v = 0.04. 


In summary, we have studied the percolation phase 
transitions of evolving random networks under a gener- 
alized Achlioptas process that the edge with minimum 
product of two connecting cluster sizes is taken with a 
probability p from two randomly chosen edges. This 
model becomes the ER network at p = 0.5_and 
the random network under the Achlioptas process at 
p= 1. Using the finite-size scaling, the percolation phase 
transitions of this model are studied. From the finite-size 
scaling forms of the largest cluster size s1(r, N; p) and the 
second largest cluster size s2(r, N; p), the fixed point of 
82/8, and the straight line of ln sı versus InN can be 
used to determine the critical points. It has been found 
that the critical points re determined from s2/s; and 
In sı are consistent and increase with p. From the slopes 
of Ins, and In(s2/s,)’ at re, we can obtain the critical 
exponent ratio 8/v and 1/v respectively. 6/v decreases 
from 0.33(1) at p = 0.5 to 0.04(1) at p = 1. 1/v increases 
from 0.33(1) at p = 0.5 to 0.50(1) p = 1. With the scaling 
variable tN‘/”, the Monte Carlo data of s2/s1(r, N; p) at 
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Figure 5: Finite-size scaling function U(tN‘/”; p) at (A) p = 
0.5 with 1/v = 0.33 and (B) p = 1.0 with 1/v = 0.5. 
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Figure 6: U (0; p) with error bars smaller than the symbol. 


a given p and different N collapse into a scaling function 


U(tN‘/";p). The percolation phase transitions of ran- 
dom networks under the GAP are continuous always and 
its universality class depends on the probability parame- 
ter p. 
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